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Abstract 

Recent progress in extending finite temperature lattice QCD 
simulations from the T axis into the finite density plane is reviewed. 
The covered topics are a calculation of the transition line by multi- 
\ dimensional reweighting, screening mass calculations in dimension- 

' ally reduced QCD and simulations at imaginary /i. 

^ ! 

^ ! 1 Introduction 

QCD at finite baryon density plays a role in nature in two rather different 
regimes: It occurs in heavy ion collisions whose initial state has non-zero 
r-| ! baryon number, and whose subsequent plasma state is of high temperature 

and low density. It is also believed to constitute the core of neutron stars, 
consisting of cold and very dense matter. These two situations correspond 
to the regions close to the axes in the tentative QCD phase diagram Fig. |[ 
Here the interest is in the former regime, whose understanding is partic- 
ularly pressing in view of experiments at SPS, LHC (CERN) and RHIC 
(BNL). 

In order to come to first principles predictions, non-perturbative meth- 
ods are required. Unfortunately, lattice QCD has so far failed to be a 
viable tool for the analysis because of the so-called "sign-problem" . The 
QCD partition function is given by 



Z = J DU det M{ij)e 



-SymIU] 
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where Sym is the pure gauge action and det M{fi) denotes the fermion de- 
terminant. For the gauge group SU(3) and /i 7^ 0, det M is complex and 
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Figure 1: Qualitative phase diagram for QCD with two massless flavours. On 
the T axis and along the dashed line there is a crossover, solid lines mark first 
order transitions. The tricritical point is expected to move closer to the T axis 



prohibits standard Monte Carlo importance sampling, for which a mani- 
festly positive measure in the functional integral is required. While still 
no ah initio solution of the sign problem is available, recent results sug- 
gest that at finite temperatures around the deconfinement transition and 
for low densities the sign problem is milder and simulations are possible. 
This contribution reviews some recent progress obtained by reweighting 
methods and simulations of dimensionally reduced QCD. It also discusses 
what we can learn from QCD with imaginary chemical potential where the 
integration measure is real and no reweighting is necessary. For a more 
complete coverage and references to other approaches see e.g. 0. 

2 Reweighting methods 

The Glasgow method [Q] evaluates the partition function by absorbing 
the complex determinant into the observable, and doing the importance 
sampling with the positive part of the measure. 



However, numerical results give for the onset of baryon density on the 
/i axis the unphysical value /xq = ?^7r/2 instead of the expected /io = 
m^/S, which is also the pathological behaviour of the quenched theory 
[Q. The reason for this failure is the poor overlap between the Monte 



for finite m^. From [Q]. 
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Carlo ensemble of configurations at /i = used to compute averages by 
reweighting, and the full ensemble at yU 7^ 0. 

This situation can be improved by splitting the determinant in mod- 
ulus and phase, det M = \ det M\ exp(i0), and incorporating the modulus 
into the measure, while only the phase is used for reweighting 



The expectation value of the phase factor can be viewed as a ratio of two 
partition functions, one of the full theory and one with only the positive 
modulus of the determinant in its measure. It thus behaves as exp(— AF), 
where the exponent is the free energy difference and an extensive quan- 
tity. Consequently, (exp(i0)) ~ exp{—V), whose measurement requires 
exponentially large statistics for realistic volumes. 

Recently considerable progress has been reported by means of a two- 
dimensional reweighting method , which in addition to /i also reweights 
in the lattice gauge coupling /?, rewriting the partition function as 



The simulation performed in for 2+1 flavours of staggered quarks is 
concerned with the critical line and the endpoint of the deconfinement 
phase transition. The simulation parameters are tuned such as to stay on 
the (pseudo-) critical line deflned by the location of Lee Yang zeros. While 
reweighting only in /i would inevitably employ an ensemble away from 
criticality, the second parameter can be used to keep the Monte Carlo 
ensemble fluctuating between the phases, as the full ensemble certainly 
would. This simulation constitutes the flrst numerical prediction of the 
(/i, T) phase diagram and is displayed in Fig. |^. However, it still has some 
caveats, and more work is required to conflrm these results. The lattices 
are still rather coarse (a ~ 0.3/m) and the volumes small. Since the sign 
problem becomes exponentially stronger with the volume, it remains to 
be seen whether inflnite volume and continuum limits can be taken with 
this method. More importantly, while the two-dimensional reweighting on 
the critical line does sample both phases, there is no guarantee that it has 
a good overlap with the full ensemble, which is still at a different point 
of parameter space. Hence it seems important to cross check these results 
by a different method. 




(3) 
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Figure 2: Phase diagram for 2+1 flavour QCD {m.,^ is about four times its 
physical value, lOniu). From [^. 



3 Imaginary chemical potential 

There are a number of suggestions to consider imaginary chemical po- 
tential, for which the integration measure is positive and simulations are 
straightforward. The connection to real chemical potential is provided by 
the canonical partition function at fixed quark number Q, which is related 
to the grand canonical partition function at fixed imaginary chemical po- 
tential by JT], § 



Z{T,Q) 



l-nT 2^ 



Z(T,fi = ifii) e 



(5) 



One proposal is to simulate Z for various values of /i/, and then numeri- 
cally do the Fourier integration 0. This becomes more and more difficult, 
however, for large Q, and extrapolation to the thermodynamic limit seems 
questionable. The method has been tested in the two-dimensional Hub- 
bard model but not for QCD. 

Another proposal, to be further pursued in the following, is to use ana- 
lyticity of the partition function to continue expectation values computed 
with Z(i/i,T). Not too far from the temperature axis (and away from a 
phase transition), physical quantities are analytic functions of fi/T having 
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a Taylor series. One can then fit the coefficients of the truncated series 
to results obtained at imaginary fi and analjd;ically continue the series to 
real /i. This has been explored for the chiral condensate in the strong 
coupling limit [0 . In the following it will be tested non-perturbatively in 
the deconfined phase in the framework of an effective theory. 



4 Dimensionally reduced QCD 

For temperatures larger than a few times the deconfinement temperature 
Tc, the static (equilibrium) physics of QCD can be accurately described by 
an effective theory, which indeed permits simulations with non-vanishing 
real quark chemical potential /i/T< 4.0, where the sign problem sets in 
1^. It permits to study any number of quarks with small or zero mass. 



The domain of parameter space where the approach is applicable contains 
the phenomenologically relevant region in which heavy ion collisions are 
operating. At and above SPS energies, densities in heavy ion collisions are 



estimated to be ixb/T <A.O |jT2|, i.e. a quark chemical potential fi/T < 1.3, 
which is well within the range where simulations are feasible. 

Equilibrium physics is described by euclidean time averages of gauge 
invariant operators, ^(x) = T Jq^^ dr A{'x, — ir). Their spatial correlation 
functions 

C(|x|) = (A(x)A(0)),~e~^-^W, (6) 

decay exponentially with distance. The "screening masses" M correspond 
to the inverse length scale over which the equilibrated medium is sensitive 
to the insertion of a static source carrying the quantum numbers of A. 

For length scales larger than the inverse temperature, |x| ~ l/gT ^ 
1/T, the time integration range becomes very small and the problem ef- 
fectively three-dimensional. The calculation of the correlation function 
C(|x|) can then be factorized: the time averaging may be performed per- 
turbatively by expanding in powers of the ratio of scales gT/T ~ g, which 
amounts to integrating out all modes with momenta ~ T and larger, 
i.e. all non-zero Matsubara modes and in particular the fermions. This 
procedure is known as dimensional reduction It is in the spirit of a 



Wilsonian renormalization group approach, where an effective action for 
coarser scales is derived by averaging over the smaller scales. The remain- 
ing correlation function of 3d fields is then to be evaluated with a 3d purely 
bosonic effective action, which describes the physics of the modes ~ gT 
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and softer. Without fermions and one dimension less, accurate infinite 
volume and continuum extrapolations of simulations are feasible. 

The effective theory emerging from hot QCD by dimensional reduction 
is the SU(3) adjoint Higgs model 0, [l^, |15| with the action 



S 



(fx 



Tr(AAo)' + mi^AK) + A3(Tr(A^ 



2\2 



(7) 



As 4d euclidean time has been integrated over, now appears as a scalar 
in the adjoint representation. The associated Hamiltonian respects SO (2) 
planar rotations, two-dimensional parity P, charge conjugation C and Aq- 
reflections i?, and screening masses can be classified by ■ 

The parameters {(/I, ms, A3} of the effective theory are via dimensional 
reduction functions of the 4d gauge coupling , the number of colours N 
and flavours A^/, the fermion masses and the temperature T. Note that 
ma ~ gT is just the leading order Debye mass. In all of the following 
fermion masses are assumed to be zero, but in principle any other values 
may be considered as well. 

With the simulation part being rather accurate, the main error in 
the correlation functions is due to the reduction step, which has been 
performed to two loops |T^. Fig. ^ compares the results for hot pure 



gauge theory as obtained in the 4d lattice theory [16| with those from the 
effective theory |[T^, |TT| . Note that the effective theory is only valid up to 
its cut-off M/T ~ Itx and above this level disagreement is to be expected. 
For the case of SU(3) one finds again quantitative agreement in the largest 
correlation length but about 20% deviation in the next shorter one. Good 
agreement between the 4d and the 3d effective theory is also found with 
other observables, like the static potential the spatial string tension 
1^ , the Debye mass |]l9l and gauge- fixed propagators pO|. Dimensional 



reduction also works to high precision from 3d finite T to 2d Thus, 
dimensionally reduced pure gauge theory gives a reasonable description of 
the largest correlation lengths in the system down to temperatures as low 



as 



2T,. 



4.1 Finite density 

The main advantage of the effective theory is that fermions, having always 
non-zero Matsubara modes, are treated analytically. Any change in the 
number of fermion species or their masses is encoded in the parameters of 
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Figure 3: Comparison of screening masses (C = +) in hot SU(2) pure gauge 
theory as determined in 4d (empty symbols) and 3d (full symbols) effective 
theory simulations. From ||l7[| . Numbers for SU(3) from [pdl |l^. 



the effective theory and only shifts the values for the screening masses in 
Fig. ^. At ~ 2 — 3Tc, the fermionic modes begin to feel the chiral phase 
transition and thus become very light through non-perturbative effects. 
At some point this effect will be so large that they constitute the lightest 
degrees of freedom and may no longer be integrated out, as demonstrated 
in a 4d simulation with Nf = 4 light fermions . 

When a chemical potential is switched on, its leading effect is to gen- 
erate one extra term in the action Eq. (|^) and to change the Debye mass. 



S 



S + iz I d^xTiAl 



T37r 



2 ' 



1 + 



3N 



f 



Wt^ 2N + N 



The extra term Sz is odd under R, C, and hence the action no longer re- 
spects these symmetries, while parity is left intact. Consequently, screen- 
ing states now are labelled by J^. 

The effective action is complex and still has a sign problem. Expecta- 
tion values have to be computed by reweighting with the complex piece 
of the action 

(Oe--^^)o 



{O) 



{cos{zSz))o ' 



(9) 
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Figure 4: Results for A'^^ = 2 and T = 2Tc. Left: Distribution of the reweight- 
ing factor (zSz) for fixed L. Right: Screening masses in the channel 0"*" T = 2Tc. 
'Scalar" states (Tr^Q, empty symbols) depend on fj./T, while "gluonic" states 
(TrF^2) full symbols ) are practically independent of it. From |11]. 



where the subscript denotes averaging with ^ = 0. Cancelling contribu- 
tions to the expectation value occur whenever (zS^) ^ 1. Fig. ^ shows 
some distributions of (zSz) as obtained by Monte Carlo for chemical poten- 
tials. As long as /i < 4T, the distribution is well contained within [— vr, tt] 
for volumes large enough so that the masses attain their infinite volume 



values without any cancellations W^ . For even larger values of the chem- 



ical potential the sign problem sets in and statistical errors explode. 

The results for the lowest screening masses are also displayed in Fig. |[ 
The different qualitative behaviour of 3d gluonic and scalar states observed 
generally in scalar gauge models p3|, 0, [lT| leads to a level crossing at 
/i/T ~ 4.0 and hence to a change in the nature of the ground state excita- 
tion. Hence there is interesting structure in the phase diagram even away 
from a phase transition, implying that the longest correlation length in 
the thermal system does not get arbitrarily short with increasing density, 
but stays at a constant level. 



4.2 Imaginary vs. real chemical potential 

With an effective theory permitting simulation of real and imaginary 
chemical potential at hand, one may now return to the suggestion of Sec. ^ 
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Figure 5: The lowest screening masses at imaginary and real /x, T = 2Tc. From 



and study the feasibiUty of analytic continuation of observables . Away 
from phase transitions, the screening masses are analytic in /i/T, as there 
are no massless modes in the theory. Moreover, since a change fi —fi 
can be compensated for by a field redefinition Aq —Aq in Eq. (|^), all 
physical observables must be even under this operation. In the 4d theory 
the same statement follows from compensating fi — > — /i by a C (or CP) 
operation. For small fJ^/T, the screening masses may thus be written as 



M 

y = Co + Ci 



nTj 



+ C2 



TlTj 



+ 



(- 

VttT 



(10) 



In the range where this ansatz fits the data of real and imaginary /i, the 
possibility of analytic continuation is easily checked by examining whether 
the Ci agree between the two data sets. 

Fig. ^ shows the lowest states in the system for real and imaginary 
potential. Note that in the imaginary /i case, data are only available 
up to ~ 1.5, and hence C2 is not well constrained. The reason is 

that imaginary chemical potentials favour a Z^r-broken minimum over the 
symmetric one in the 4d effective Aq potential, once |/i|/T > 7r/3. Thus a 
phase transition occurs and analyticity is lost 0. 

In the region up to this critical value, however, one observes firstly that 
the coefficients Co,i are sufficient to fit the data, and secondly that they are 
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completely compatible between the real and imaginary fi data sets [p^ . 
Having only two coefficients is equivalent to susceptibility measurements 
at /X = 0. The ffist conclusion then endorses work based on that latter 
approach ||2^. The second conclusion encourages us to consider 4d QCD 
at imaginary /i and study the vicinity of the critical line. 



5 QCD at imaginary chemical potential 

From the fact that screening masses decrease with imaginary /i it follows 
that Tc(/m(/i)) is an increasing function, leading to Fig. ^ (left) as the 
complex generalization of the phase diagram Fig. First results of 




5.40 



Figure 6: Left: Qualitative phase diagram for complex /i. Right: Susceptibility 
of the Polyakov loop for imaginary fj,. T increases with (3. 



simulating two-flavour QCD with staggered fermions at imaginary /i [26 



fully conffim this picture, as shown in Fig. ^ (right). The peak of the 
susceptibility of a gauge invariant operator occurs at a (pseudo-) critical 
point. The figure clearly shows that the critical coupling (3c (and with it 
Tc) is growing with Furthermore, on finite volumes susceptibilities are 
analytic funtions. This means that Tc(/i) itself is analytic and it might be 
possible to continue it from imaginary to real /i. This question is currently 
under investigation. 
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6 Conclusions 



Recent developments have shown that at high temperature, despite the 
sign problem, lattice simulations can be extended away from the tempera- 
ture axis into the finite density plane. Considering dimensionally reduced 
QCD gives a detailed picture of screening masses in the deconfined phase, 
covering the experimentally relevant parameter space. Furthermore, it 
establishes that observables depend only weakly on /i and hence suscep- 
tibility measurements as well as analytic continuation from imaginary to 
real /i are possible. Multidimensional reweighting has even produced a 
result for the critical line of the deconfinement transition and its end- 
point. Because it uses reweighting, however, an independent check would 
be desirable. 
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